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Abstract. 

We present a probabilistic approach to water- water hydrogen bonding that allows one to obtain an analytic 
expression for the number of bonds per water molecule as a function of both its distance to a hydrophobic 
particle and hydrophobe radius. This approach can be used in the density functional theory (DFT) and 
computer simulations to examine particle size effects on the hydration of particles and on their solvent-mediated 
interaction. For example, it allows one to explicitly identify a water hydrogen bond contribution to the external 
potential whereto a water molecule is subjected near a hydrophobe. The DFT implementation of the model 
predicts the hydration free energy per unit area of a spherical hydrophobe to be sharply sensitive to the 
hydropobe radius for small radii and weakly sensitive thereto for large ones; this corroborates the vision of 
the hydration of small and large length-scale particles as occurring via different mechanisms. On the other 
hand, the model predicts that the hydration of even apolar particles of small enough radii may become 
thermodynamically favorable owing to the interplay of the energies of pairwise (dispersion) water-water and 
water-hydrophobe interactions. This sheds light on previous counterintuitive observations (both theoretical 
and simulational) that two inert gas molecules would prefer to form a solvent-separated pair rather than a 
contact one. 
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1 Introduction 

A particle whereof the accommodation in water is accompanied by an increase in an associated free 
energy is called "hydrophobic" and is often referred to as a hydrophobe. The thermodynamically 
unfavorable dissolution of a hydrophobe (whether microscopic or macroscopic) is hydrophobic hy- 
dration; the corresponding free increase results from structural (and possibly energetic) changes in 
water around the hydrophobe. The total volume of water affected by two hydrophobes is smaller 
when they are close together than when they far away from each other. This gives rise to an effective, 
solvent-mediated attraction between them which is also referred to as hydrophobic attraction. 

Hydrophobic effects (hydration and attraction) play a crucial role in various physical, chemical, 
and biological phenomena. 1 " 4 They also play an important role in the formation, stability, and 
unfolding of the native structure of a biologically active protein which constitute the core of two 
most exciting (and intrinsically related) topics of modern biophysics, namely, "protein folding" and 
"protein denaturation" , 5 ' 6 although an assortment of interactions (including those of hydrophilic 
character 7 ) would most likely determine the driving force of these amazing phenomena. 

Various mechanisms have been suggested to understand hydrophobic effects at a fundamental 
level and develop a general theory of hydrophobicity. 8-11 Virtually all theoretical models involve the 
hydrogen bonding ability of water as a key element. 

The structure of liquid water, its dependence on the external conditions, and the role of structural 
changes in hydrophobic phenomena have long been the subject of intense research. For ambient 
conditions, it was first described as a locally ordered tetrahedral network of water molecules. 12 Various 
anomalous properties of water (such as the density maximum at 4° C at atmospheric pressure, local 
maximum and minimum of the isobaric heat capacity at constant pressure, etc..) are attributed to 
the ability of its molecules to form hydrogen bonds, strong directional bonds with energy much larger 
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than the thermal energy ksT. As an example of water structure effects on hydrophobic phenomena, 
the biological activity of proteins appears to depend on the formation, maintenance, and breakup of 
a 2-D hydrogen-bonded network spanning most of the protein surface and connecting all the surface 
hydrogen-bonded water clusters. 13 

Although much experimental, computational, and theoretical research has been carried out, many 
thermodynamic and molecular aspects of hydration remain to be clarified. In the positive (unfavor- 
able) free energy of dissolving a hydrophobe, the positive entropic contribution (due to the negative 
entropy change) dominates over the enthalpic contribution at room temperatures. The total con- 
tribution of hydrogen bonding to the hydration enthalpy depends both on the single bond energy 
and the number of bonds that a water molecule can form in the hydrophobe vicinity and in the 
bulk. The propensity of water to form all possible hydrogen bonds, on one hand, and the constraint, 
imposed by a hydrophobe on the configurational space available to vicinal water molecules, on the 
other hand, lead to a large entropic cost. Some investigations have suggested that water is more 
structured near a hydrophobe, with water-water hydrogen bonds both labile and stronger than in 
bulk, but several experimental and theoretical studies have reported the opposite results. Despite 
remaining controversies, the dependence of hydrophobic phenomena on the length scales of solute 
particles is considered to be proven. 14-16 

The hydration of small hydrophobic molecules (of sizes comparable to a water molecule) is believed 
to be entropically "driven" (and so is their solvent-mediated interaction). 10 ' 11 Such molecules can fit 
into the water hydrogen-bond network without destroying any bonds. While this results in a negligible 
enthalpy of hydration, the solute constrains some degrees of freedom of neighboring water molecules 
which gives rise to negative hydration entropy and hence to positive hydration free energy. However, 
such a simple mechanism has recently come under scrutiny 10 ' 11 ' 14 ' 15 because there are simulations 17 ' 18 
and theory 19 suggesting that, under some conditions, the hydration of small hydrophobic molecules 
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could be entropically favorable. 

The hydration of large hydrophobic particles is believed to occur via a different mechanism. 10,11 ' 19,2( 
When inserted into liquid water, a large hydrophobe breaks some hydrogen bonds in its vicinity. This 
would result in large positive hydration enthalpy and hence in a free energy change proportional to 
the solute surface area (as opposed to being proportional to the solute volume for small hydrophobes). 
Thus, the hydration of large hydrophobic particles is expected to be enthalpically driven (and so is 
their solvent-mediated interaction). 

As the thermodynamics of hydration is expected to change gradually from entropic for small 
solutes to enthalpic for large solutes, so are the structural properties of liquid water in the vicinity of 
the solutes. It was argued 21 that if the solute-water attraction is sufficiently weak, there may exist 
a thin film of water vapor near large hydrophobic solutes but not small ones. This generated much 
controversy. 10 ' 11 ' 18 ' 20 ' 22 

Hereafter we present a model for water-water hydrogen bonding that allows one to obtain an 
analytic expression for the number of bonds per water molecule as a function of both its distance to 
a hydrophobe and hydrophobe radius. This function can be used in the density functional theory 
(DFT) and computer simulations (either Monte Carlo or Molecular Dynamics) to examine particle 
size effects on the hydration of particles and on their solvent-mediated interactions over the entire 
small-to-large length-scale range. 

Note that we do not investigate drying or wetting transitions 23 as such; once the accomodation of 
the solute particle in liquid water occurred, the fluid density distribution in the vicinity of the solute 
and in the entire system is not subject to any "transformations". The latter may be induced only by 
changes in external thermodynamic variables (either temperature or pressure or chemical potential). 
We will consider a "static" version of the hydration phenomenon, wherein the state of the system 
does not change after hydration occurred. 
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2 The number of hydrogen bonds per water molecule near a spher- 
ical hydrophobic surface 

Consider a spherical hydrophobic particle of radius R immersed in liquid water (Figure 1). Even if 
one assumes that the intrinsic hydrogen bonding ability of a water molecule is not affected by the 
hydrophobe, in its vicinity a "boundary" water molecule forms a smaller number of bonds than in bulk 
because the surface restricts the configurational space available to other water molecules necessary 
for a boundary water molecule to form hydrogen bonds. The probabilistic model allows one to obtain 
an analytic expression for the average number of bonds that a boundary water molecule can form as a 
function of its distance to the hydrophobe and hydrophobe radius. A boundary hydrogen bond may 
be slightly altered energetically compared to the bulk one, but such alteration is still uncertain 24 " 26 
and will be neglected hereafter. 

In the probabilistic hydrogen bond (PHB) approach, 27 a water molecule is considered to have four 
arms each capable of forming a single hydrogen bond. The configuration of four hydrogen-bonding 
(hb) arms is rigid and symmetric (tetrahedral) with the inter-arm angles a = 109.47° (Fig.l). Each 
hb-arm can adopt a continuum of orientations subject to the constraint of tetrahedral rigidity. A 
water molecule can form a hydrogen bond with another molecule only when the tip of any of its 
hb-arms coincides with the second molecule. The length of a hb-arm thus equals the length of a 
hydrogen bond rj, assumed independent of whether the molecules are in bulk or near a hydrophobe. 
The characteristic length of pairwise interactions between water molecules and molecules constituting 
the hydrophobe is also assumed to be rj. 

The location of a water molecule is determined by the distance r from its center to the center of 
the hydrophobe which is also chosen as the origin of the spherical coordinate system. The distance 
x between water molecule and hydrophobe is defined as x = r — R (Fig.l). 
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Denote the number of hydrogen bonds per bulk water molecule by n& and the average number of 
hydrogen bonds per boundary water molecule by n s . The latter is a function of radius R and distance 
X . 1.6. j — Tl § (R,x). If x > 2r], the number of hydrogen bonds that the water molecule can form 
is assumed to be unaffected by the hydrophobe: n s (R,x) = nf, for x > 2rj. On the other hand, the 
function n s (R,x) attains its minimum at x = 77, because at this distance the configurational space 
available for neighboring water molecules is most restricted compared to bulk water. A spherical 
layer of thickness r\ from r = R + rj to r = R + 2n is referred to as the solute hydration layer (SHL). 

In the spirit of the PHB approach 27 let us represent the function n s = n s (R,x) as 

n s = kibi + k 2 b\ + k 3 bf + k 4 bf, (1) 

where b\ is the probability that one of the hb-arms (of a bulk water molecule) can form a hydrogen 
bond and the coefficients k\, k 2 , ks, and k 4 depend on R and x, and so does n s . Equation (1) assumes 
that the intrinsic hydrogen-bonding ability of a water molecule (the tetrahedral configuration of its 
hb-arms and their lengths and energies) is unaffected by the hydrophobe. 

The functions k\ = k\(R, x), k 2 = k 2 (R,x),k3 = ks(R,x), and k 4 = k 4 (R,x) can be evaluated by 
using geometric considerations (see the Appendix). They all become equal to 1 at x > 2rj, where 
eq.(l) reduces to its bulk analog, rib = b\ + b\ + b\ + b\ (see the Appendix). Since experimental data 
on nb are readily available, one can find b\ as a positive solution (satisfying < b\ < 1) of the latter 
equation. 

Thus, equation (1) provides an efficient pathway to n s as a function of x and R. It takes into 
account the constraint that near the hydrophobe some orientations of the hb-arms of a boundary 
water molecule cannot lead to the formation of hydrogen bonds. This constraint depends on the 
distance betwe water molecule and hydrophobe and on the hydrophobe radius, whence the R- and 
x-dependence of k\, k 2 , and k 4 . 
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Figure 2 presents the function n s (R, x) for a spherical hydrophobe immersed in water at temper- 
ature T = 293.15 K, which corresponds to n& = 3.65 hence b\ = 0.963707. In Fig. 2a, n s is plotted 
vs £ = (x/77 — 1) for various radii R. As expected,n s monotonically increases from its minimum at 
x = i] to its maximum bulk value n& at x = 2r/. For a flat hydrophobic surface (i? = 00), molecular 
dynamics simulations 28 ' 29 previously reported such behavior of n s (although with some oscillations in 
ref.29). In Fig. 2b, n s is shown as a function of R at different distances x. At any x, n s monotonically 
decreases from its maximum for the smallest particle R = to its minimum for the largest particle 
(i? = 00). Besides, for any fixed x, as R increases from to 00, n s approaches its asymptotic value 
for a flat hydrophobic surface, n s (oo,x), for particles of radii as small as R « 30??. 

3 Implementation of the probabilistic hydrogen bond model in the 
density functional theory 

The fluid density distribution near a rigid surface can be efficiently studied by using computer 
simulations or DFT. 30-32 As an illustration of the PHB approach, let us implement it into DFT. 
The latter usually treats the interaction of fluid molecules with a foreign (impenetrable) substrate 
in the mean-field approximation whereby every fluid molecule is considered to be subjected to an 
external potential, due to its pairwise interactions with the substrate molecules. 31,32 The substrate 
effect on the ability of fluid (water) molecules to form hydrogen bonds had been previously neglected. 
However, using the PHB model, one can explicitly implement that effect in the DFT formalism and 
clarify its role in the length-scale dependence of hydrophobic hydration. 

To apply DFT to the thermodynamics of hydrophobic phenomena, it is necessary to know the total 
external potential field J7 ext = U cxt (R,x) whereto a water molecule is subjected near a hydrophobic 
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particle. This potential can be written as 

[/ cxt = [/ p xt + ce xt , (2) 

where ?7 e h xt = C/ e h xt (R, x) is the water- water hydrogen bond contribution to f/ oxt , and [7 c p xt = U^ t (R,x) 
represents the external pairwise potential exerted by all the molecules constituting the hydrophobe 
on a water molecule. 

While various models were designed 31-33 for C/ C p xt , the hydrogen bond contribution C/ e h xt had been 
conventionally neglected until recently 34 ' 35 . This contribution, C/ h xt , is due to the deviation of n s 
from rib as well as the (possible) deviation of e s from e b (the latter effect is neglected hereafter). It 
can be determined as 

Utt = \{tsn s -e b n b ). (3) 

The first term on the RHS of eq.(3) represents the total energy of hydrogen bonds of a water 
molecule at a distance x from the surface of a particle of radius R, whereas the second term is 
the energy of its hydrogen bonds in bulk (at x — > oo); the factor 1/2 is needed to prevent double 
counting the energy because every hydrogen bond and its energy, either e s or e b , are shared between 
two molecules (in refs.34 and 35 the analogous equation for a planar surface was mistyped, as the 
factor 1/2 was missing). Note that U^ t (R,x) ^ only for r] < x < 2rj. 

In DFT, the grand thermodynamic potential f2 of a nonuniform single component fluid, subjected 
to an external potential U ext (representing the hydrophobe), is a functional of the number density 
p(r) of fluid molecules 

n[p(r)] = T h [p(r)] + \j j drdv'p(vW)^ t (\v-v'\) 

+ JdrU cxt (R,r)p(r)-pjdrp(r), (4) 

where ^[^(r)] is the intrinsic Helmholtz free energy functional of hard sphere fluid, p is the chemical 
potential, and (/> at (|r— r'|) is the attractive part of the interaction potential between two fluid molecules 



located at r and r'; the integrals are taken over the volume V of the system. Among various models 
for T h [p(r)], the weighted density approximation (WDA) 30 ' 36 ' 37 with a weight function independent 
of weighted density represents an optimal combination of accuracy and simplicity. It is non-local 
with respect to p(r); it takes into account short-ranged correlations and captures the fluid density 
oscillations near a hard wall. We hereafter adopt the WDA version of DFT. 

The key element of WDA is the weighted density p(r) determined in terms of p(r) via an implicit 
equation 

p{r) = J dr' p(t')w(\t' - r|;p(r)), (5) 

where w{\r' — r|;p(r)) is the weight function. Although in more sophisticated versions 30 ' 36 of WDA 
w(\y' — r|;p(r)) depends on p(r), we will hereafter adopt its simpler version wherein the weight 
function is independent 30 ' 37 of p(r). 

For hydrophobic hydration in an open system of constant p, V, and T (grand canonical ensem- 
ble), the equilibrium density profile is obtained by minimizing Q[p(r)} with respect to p(r). The 
corresponding Euler-Lagrange equation can be written as 

p = k B Tln(A 3 p(r)) + W(r;p(r)), (6) 

where A = (h? /2itmkBT) l l 2 is the thermal de Broglie wavelength of a molecule of mass m (h and 
ks being Planck's and Boltzmann's constants) and W(r; p(r)) is a function of r and a functional of 

p(r) .28-30 

W(t-p(y)) = UUr) + jdr'p(r')U\r-r'\) 

+ AMp(r)) + jdr'p(r')A^ h (p(r'))w(\r'-r\). (7) 

Here a (|r — r'|) is the attractive part of the interaction potential between two fluid molecules located 
at r and r', whereas Aip h (p) is the configurational part of the free energy of hard sphere fluid per 
molecule, with Aip' h (p) = dAip h (p)/dp. 
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The hydrophobe being spherical, the external potential is a function of a single variable x = 
r — R, and the equilibrium density profile obtained from eq.(6) is a function of a single variable r: 
p(r) = p(r). The substitution of p(r) into eq.(4) provides the grand thermodynamic potential Q of 
the non-uniform fluid with a hydrophobe therein. The grand canonical free energy of hydration is 
AS1 = Q, — Qq, where Qo is the grand thermodynamic potential of uniform liquid water without a 
hydrophobe therein. 

4 Numerical Calculations 

For a numerical illustration, we considered the hydration of a spherical hydrophobe (taking R/rj = 
1,3,5,7,10,15,20,30,50,100) in the model water at T = 293.15 K and p = -11.5989 k B T corre- 
sponding to its two-phase equilibrium. The liquid state of bulk water was ensured by imposing the 
appropriate boundary condition onto eq.(6), p(x) — >■ p\ as x — > oo, with p\ the bulk liquid density. 
The densities p v and p\ of coexisting vapor and liquid, respectively, are determined by solving the 
equations p(p,T)\ p=f)v = p(p,T)\ p=pi , p(p,T)\ p=f!v = p(p,T)\ p=pi , requiring the chemical 

potential p = p(p, T) and pressure p = p(p, T) to be the same throughout both coexisting phases. 

The chemical potential of a uniform hard sphere fluid p h and the configurational part Aip h = 
Aijj h (p, T) of the free energy of a hard sphere fluid were modeled in the Carnahan-Starling approximation, 31 ' 32 ' 
whereas for the weight function w(\r' — r|; p(v)) in eqs.(5),(7) we adopted a /5-independent version 37 

£U — 3£ ) 3 
(1 - £y vrr? 4 

with 6 (it) being the Heaviside (unit-step) function. 

The pairwise interactions of water molecules were modeled by using the Lennard-Jones (LJ) po- 
tential with the energy parameter e ww = 3.79 x 10~ 14 erg and the diameter d of a model molecule 
set to be r\. The attractive part at of pairwise water-water interactions was modeled via the Weeks- 
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Chandler-Anderson perturbation scheme. The interaction potential between water molecule and 
molecule of a hydrophobe was assumed to be of LJ type with an energy parameter e wp and a length 
parameter 77. Integrating this interaction with respect to the position of the molecule of the hy- 
drophobe over the hydrophobe volume Vr = 47ri? 3 /3, one can obtain the pairwise contribution C/ e p xt 
into f7 eX f We assumed the dimensionless number density of molecules in the hydrophobe to be 
p p rf ss 1 and considered five values for e wp /e ww (0.6, 0.66, 0.75, 0.9, 1.2) to mimic various degrees 
of hydrophobicity. The density profiles and free energies of hydration thus obtained are shown in 
Figures 3-5. 

Figure 3 presents the density profiles near a spherical hydrophobe of radius R for R/r] = 
1,3,5,7,10,15,30,100 and e wp /e ww = 0.75 . The profiles in this Figure (as well as in Figure 4) 
are presented in a coordinate system with the origin r = located in the center of the hydrophobe; 
there are no water molecules at r < R+r] (the space < r < R is occupied by the hydrophobic sphere 
and the layer R <r < R + r] is excluded to fluid molecules). As clear, the hydrophobe radius greatly 
affects the distribution of vicinal water molecules. The oscillations in the density profile gradually 
disappear as R increases. They are well pronounced for R/r] = 1, but virtually non-existent for 
particles R/r] > 7. As R increases, a thin depletion layer around the particle (virtually non-existent 
for R/r] = 1) becomes more developed, with its density approaching that of vapor and its thickness 
approaching 77. 

This is consistent with the largely accepted wisdom concerning the much discussed issue whether 
or not there is a vapor-like layer near a large hydrophobe in liquid water. As now widely agreed 
upon, even if (and when) such a layer exists, it should be expected to be of molecular thickness 
only. 10 ' 11 ' 19 ' 21 " 23 

Furthermore, the behavior of fluid density profiles in Fig.3 is consistent with our previous finding 34 '' 
that the hydrogen bond contribution to the external potential plays a crucial role in the formation 
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of a thin "strong depletion" layer (of density much lower than liquid and of thickness of a molecular 
diameter) between liquid water and planar hydrophobic surface even for weakly hydrophobic sur- 
faces (with high e wp /e ww ). Indeed, as R increases, the geometric constraint on the ability of a vicinal 
water molecule to form hydrogen bonds strengthens, the repulsive contribution U^ xt to c7 oxt increases, 
whence the widening and stronger depletion of the vicinal water layer around the hydrophobe if the 
"pairwise" hydrophobe- water (attractive) contribution [/f xt (determined by the ratio e wp /e ww ) is not 
too large (by absolute value). 

In order to further clarify the effects of water-hydrophobe attraction and hydrophobe radius on 
fluid (water) density profiles, they are plotted in Figure 4 for three different values of e wp /e ww and 
two different radii R, namely, R/rj = 1 (Fig.4a) and R/rj = 15 (Fig.4b). Three profiles shown in 
each of Figs.4a and 4b correspond to e wp /e ww = 0.6,0.75,1.2 (from bottom to top, respectively). 
As clear, the strengthening of pairwise intermolecular fluid-hydrophobe interactions by 100% has a 
little effect on density profiles near a sufficiently large hydrophobe {R/rj = 15); the thickness of the 
depletion layer remains virtually unaffected (roughly equal to rj) and the fluid density therein remains 
several orders of magnitude lower than its bulk value (i.e., the depletion layer remains vapor-like). 
On the other hand, for a molecular size hydrophobe (R/n = 1), the increase of e wp /e ww by 100% 
leads to a drastic change in the nature of the water depletion layer near the hydrophobe; it becomes 
significantly narrower and from being a vapor-like one transforms into a liquid-like one. These results 
are in qualitative agreement with the previously reported ones obtained via molecular dynamics 
simulations 40 of the SPCE water model and via Monte Carlo simulations 41 of the TIP4P water 
model. The latter study also reported the analogous behavior for all fluids, including nonassociating 
ones (without hydrogen-bonding ability). That is somewhat dissimilar from our previous studies 34 
showing that that even for a relatively strong hydrophobic planar surface (with low e wp /e ww ) the 
conventional contribution to the external potential (due to pairwise interactions between a water 
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molecule and those of the substrate) cannot cause the formation of a vapor-like layer near the surface, 
although it does lead to the formation of a depletion layer with a weak decrease in the vicinal fluid 
density compared to the bulk one. 

In the framework of the proposed approach, the thickness of and the density in the depletion layer 
are determined by the interplay of two effects. On one hand, the constraint on the configurational 
space available to water molecules, wherewith a selected molecule can form hydrogen bonds, is 
taken account of via the function n s (R,x), eq.(l). On the other hand, the model incorporates 
attractive interactions between water molecule and hydrophobe, eq.(2), whereof the strength can be 
characterised by the positive parameter e wp . For a given thermodynamic state of the system and 
the nature of the hydrophobe (represented by £ wp ), the result of this interplay naturally depends on 
the hydrophobe size (radius). When the former effect predominates over the latter, the depletion 
layer is vapor-like. Otherwise, for relatively weakly hydrophobic particles (that are not too large), 
the vicinal water layer is just slightly depleted compared to the bulk liquid. 

Figure 5a presents the grand canonical free energy of hydrophobic hydration AO as a function 
of the hydrophobe radius R (note that the curves are provided only for guiding the eye; the actual 
calculated points are at R/r) = 1,3,5,7,10,15,20,30,50, and 100). The intrinsic hydrophobicity of 
the particles is assumed to be independent of R, with e wp /e ww = 0.75. The hydration free energy is 
expressed in units of ksT per "dimensionless unit area"; the dimensionless AO in Figure 4 is obtained 
by dividing AO by ksT and by 4ttR 2 /rj 2 . The variable sensitivity of AO to R is a clear indication 
that the hydration of small and large length-scale particles occurs via different mechanisms and that 
the hydrogen bond contribution to U cxt plays a key role in this process. The model predictions for 
AO for small i?'s are consistent with the experimental data on the hydration free energy of methane, 
ethane, propane, and n-butane at the temperature T = 300 K, as compiled in ref.8; considering a 
methane molecule as a sphere and ethane, propane, and n-butane molecules as cylinders, one can 
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roughly estimate the experimental A£l to be 0.6 for methane, 0.4 for ethane, and 0.3 for propane 
and n-butane. 

The dominant role of the hydrogen bond network in hydrophobic hydration is emphasized in 
Figure 5b where Af2 is plotted vs e wp /e ww for various radii R. Each curve in Fig.5b corresponds 
to a fixed R, with R/rj = 1,3,7,15,30 from bottom to top (again the curves are provided only 
for guiding the eye; the actual calculated points are at e wp /e ww = 0.6,0.66,0.75,0.9, and 1.2). As 
expected, the hydration free energy per unit area decreases with increasing degree of hydrophobicity 
and, for small enough particles, Ail may even become negative. The model predictions suggest 
that the hydration of even apolar particles of radii R < 3r? may be thermodynamically favorable 
(the hydration free energy being negative) if the pairwise (LJ-type) interactions between a water 
molecule and a molecule constituting the hydrophobe are comparable with or stronger than the 
(LJ-type) interactions between two water molecules, i.e., if e wp /e ww > 1. This result clarifies some 
previous simulational and theoretical observations 18-20 that two inert gas molecules would prefer to 
form a solvent-separated pair rather than a contact pair (dimer). 

5 Conclusions 

Concluding, we emphasize that the PHB approach allows one to obtain an analytic expression for 
the average number of hydrogen bonds per water molecule near a spherical hydrophobic particle 
as a function of both the particle radius and the distance between water molecule and hydrophobe 
surface. This function can serve as a foundation for elucidating various aspects of hydrophobic 
phenomena, particularly their length-scale dependence, either via computer simulations (Monte Carlo 
and Molecular Dynamics) or DFT. For example, this function allows one to explicitly identify an 
additional contribution to the external potential exerted by the hydrophobe on a water molecule; 
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it is due to the alteration of water hydrogen bonding near a hydrophobe. Thus, one can efficiently 
implement the hydrogen bonding ability of water molecules in DFT to examine the particle size 
dependence of hydrophobic hydration. 

As a numerical illustration of the combined PHB/DFT approach, we have studied the hydration 
of spherical particles of various radii and various hydrophobicity in a model water. The numerical 
results for the hydration free energy of small size hydrophobes are consistent with the experimental 
data on the hydration of small alkanes. The free energy of hydration per unit area of a spherical 
particles is predicted to have a varying sensitivity to the particle radius which is a clear indication 
that the hydration of small and large length-scale particles occurs via different mechanisms. On 
the other hand, the model predictions suggest that the hydration of even apolar particles of small 
enough radii may be thermodynamically favorable if the pairwise (dispersion) attraction between 
a water molecule and a molecule constituting the hydrophobe is comparable with or stronger than 
the pairwise (dispersion) attractions between two water molecules. This result at least partially 
clarifies some previous simulational and theoretical observations (rather counterintuitive from the 
conventional point of view on hydrophobicity) that two inert gas molecules would prefer to form a 
solvent-separated pair rather than a contact pair (dimer). 

Note that in the PHB approach the tetrahedral rigidity (both geometric and energetic) of the 
hb-arms of a water molecule is assumed only for the analytical simplicity. It can be eliminated 
to allow for the geometric deformation and energetic alteration of this configuration depending on 
the sequence in which the hb-arms are engaged or due to the proximity of the water molecule to 
the hydrophobe. These modifications will just render the model more complicated for analytical 
treatment and can be expected to relatively weakly affect the model predictions. 
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Appendix. Derivation of the coefficients ki, &2, ^3, and k± as functions 
of R and x 

To find the function n s (R,x), first consider its bulk analog, n&, and represent it as (see ref.[27] in the 
main text) 

n b = bi + 62(1) + 63(2,1) + fe 4(3,2,i) . (Al) 

where 61 is the probability that one of the hb-arms (of a bulk water molecule) can form a hydrogen 
bond, 6 2 (i) is the probability that a second hb-arm can form a hydrogen bond subject to the condition 
that one of the hb-arms has already formed a bond, 63(2,1) is the probability that a third hb-arm can 
form a hydrogen bond subject to the condition that two of the hb-arms have already formed bonds, 
and 64(3,2,1) is the probability that the fourth hb-arm can form a bond subject to the condition that 
three of the hb-arms have already formed bonds. 

Note that the probability 61 can be formally represented as a product b\ = Ps->nPn^>s, where 
Ps-¥N is the probability that the tip of any hb-arm of molecule S roughly coincides with molecule 
N and Pn^s is the probability that the tip of any hb-arm of molecule N roughly coincides with 
molecule S. (Similar considerations are valid for 6 2 (i), 6 3 ( 21 ), and 64(3^,1) as well). Neither Ps-+n 
nor Ps^s can be found in the framework of our simple model, but their product (i.e., 61) can be 
determined from readily available experimental and simulational data on n^. 

Indeed, in the chosen model of a water molecule the events of formation of bonds by the hb- 
arms (in bulk water) can be considered as independent of each other, so that 6 2 (i) = 6?, &3(2,i) = 
6?) 64(3 2,1) = b\. Thus, the probability 61 can be evaluated as the positive solution of the equation 
nb = 61 + b\ + b\ + bf satisfying the condition < b\ < 1. The latter representation of n& implies 
that the intrinsic hydrogen bonding ability of each arm is independent of whether the other arms 
have been already engaged in hydrogen bonds or not. That is, when the first hb-arm of a water 
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molecule forms an actual bond, the electron density distribution in a water molecule determining 
the ability of the other three hb-arms (that are not engaged yet) to form bonds (and their potential 
orientations) remains unaffected (there is no issue with the availability of water molecules necessary 
for the selected bulk molecule to form bonds). 

For a boundary water molecule, let us represent n s in a form: 

n s = si + s 2 (i) + $3(2,1) + $4(3,2,1)- (A2) 

Here si = si(R,x), s 2 (i) = s 2 (i)(R,x), s 3(2i i) = 33(2,1) s 4 ^ 2 ,i) = $4(3,2,1) are proba- 

bilities analogous to b±, 6 2 (i), ^3(2,1), ^4(3,2,1) subject to the constraint that some orientations of the 
hb-arms cannot lead to the formation of hydrogen bonds because of the proximity to the hydropho- 
bic particle. The severity of this constraint depends on the distance of the water molecule to the 
particle, hence the x'-dependence of s±, $ 2 m, $3(2,1), $4(3,2,1)- Again, as a first approximation the 
intrinsic hydrogen-bonding ability of a water molecule (i.e., the tetrahedral configuration of its hb- 
arms and their lengths and energies) can be considered to be unaffected by its proximity to the 
hydrophobic particle so that the latter only restricts the configurational space available to other 
water molecules necessary for this boundary water molecule to form hydrogen bonds. Thus, one can 
relate s l: s 2 (i), s 3 ( 2 ,i), and s 4 (3,2,i) to h ^ b 2(i), b 3(2,i), and 64(3,2,1), respectively, as 

Si = kih, S 2 (l) = fc 2 & 2 (l), S 3(2,l) = ^3^3(2,1), s 4(3,2,l) = ^4^4(3,2,1), (A3) 

where the coefficients fci, k 2 , A; 3 , and k 4 are functions of R and x (with their dependence on the 
boundary water molecule orientations averaged) and can be evaluated by using geometric consider- 
ations. 
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Coefficient k\ 

The coefficient k\ is calculated by taking into account that a boundary water molecule can form a 
hydrogen bond "almost" like a bulk molecule except for the constraint that the tip of the hb-arm 
(arm 1) must not be too close to the surface of the hydrophobic particle of radius R. Select an 
arbitrary water molecule at a distance 77 < x < 2rf from that surface and denote it S (Figure 4). 
Any of its hb-arms can form a hydrogen bond if the tip of the arm is located anywhere on a sphere 
of radius r\ (centered at S) from which a spherical cap is cut out by the sphere ib of radius R + rj, 
inner boundary of the SHL of particle R. Denoting the corresponding solid angle by Q s (R,x), one 
can write 

f ef x (x) 

n s (R,x) = 2ir / dGi sin(Gi), (A4) 

Jo 

where 0i is the angle between hb-arm 1 and radial axis r (with the origin in the center of the 
hydrophobe and passing through the molecule S), with Qf^.(R, x) = arccos[— (2(i? + rj + x) + rj 2 + 
x 2 ) /2{R + ri + x)rj\ is the maximum angle ©1 at which hb-arm 1 can still form a bond. The probability 
si(R, x) that any one of hb-arms of molecule S can form a hydrogen bond is related to 61 via 

Sl (R,x)= ' ' V (A5) 

where = 4tt. Integrating the RHS of eq.(S4), substituting the result into eq.(S5), and taking into 
account eq.(S3), one obtains the coefficient k\ = k±(R,x) to be : 

, _ I L (2(R + y)x + V 2 + x 2 ) \ 

kl ~2{ 1+ 2(R + V + x) V J' (A6) 

Coefficient ^ 

The coefficient ki is calculated by assuming that hb-arm 1 has already formed a bond in an arbitrary 
orientation 0! (with respect to the radial axis r) and by taking into account that hb-arm 2 can form 
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a bond subject to the condition that its tip must not lie closer to the particle surface than the SHL 
inner boundary ib of hydrophobe R. 

For a bulk water molecule the probability &2(i)> that the second hydrogen bond forms once the 
first one has formed is proportional to the full length L& of the circle C*2(l) formed by the possible 
loci of the tip of the engaged second hb-arm (hb-arm 2) subject to the restriction that the angle 
between the two hb-arms remains a. Since the radius of that circle is equal to Rc 2 m = 7 ?sin(a), we 
have 

L b = 2TT7]s'm(a). (A7) 

However, for the molecule S the possible loci of the tip of engaged hb-arm 2 (subject to the 
restriction that the angle between it and hb-arm 1 is a) constitute just a part of the circle of radius 
r/sin(a), - the other part (a circular arc) being excluded by the proximity of the SHL inner boundary 
ib. The length of the "available" part of this circle is a function of R, x, and 6i and will be denoted 
by L s = L s (R,x,@i) (by definition, L S (R, x, 0i) = L5 for x > 2rj). The probability sL^ that 
molecule S engages in a second hydrogen bond once its hb-arm 1 has already formed a bond (at 
given R, x, 61) is proportional to L s : 

4(1) = 4(1) (R,x,Qi) = Ls{R, * b ,@l) b 2{1) (A8) 

Let us introduce the Cartesian coordinate system with the origin O in the center of particle R, 
axis z coinciding with the radial axis r, and axis y directed from the origin towards the projection of 
the tip of hb-arm 1 onto the x — y plane. Depending on the angle 0i, distance x, and radius R, the 
circle C2(l) may either intersect the SHL inner boundary (which is a sphere of radius R + rj centered 
at 0, hereafter denoted S l $ HL ) or not. In the former case, there may be either two intersection points 
which can degenerate into one point in the limiting at some particular orientation for given R and 
x. Assuming that there are two intersection points of sphere Sl^^ and circle C2(l), let us denote 
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their Cartesian coordinates by x_, y_, z_ and x+, y+, z+. Clearly, y± = y_ = y + , z± = = z+, and 
x_ = — x + ; one can choose the notation so that x + > 0; clearly the coordinates of both points are 
functions of R, x, ©i. 

Further, let us define the angle ©io(-R, x) to be the value of the angle 0i at which the circle C2(l) 
just "touches" the sphere S l g HL , and introduce x° = x°(R) and R° as the solutions of equations 
Q\o(R, x°) = (with respect to x°) and x°(R°) = (with respect to R°), respectively. One can thus 
obtain R° » 0.50009, 



(A10) 



x° = -R-n cos(q) + ^R 2 + 2Rr] + r/ 2 cos 2 (a), e 10 (R, x) = -a + Of x {R, x), (A9) 

For < R < R°, one can show that 

L b if 0i G [0; 010(12, x)] 

Lf if ©i € }& 10 (R,x)-&f M (R,x)} 

L b if 0i G }ef M (R,x);@f I x (R,x)] 

where 0^ M (i?,x) = min[0^( J R, x), 0^(i?, x)], &? x (R,x) = 2tt - a - @^ X {R, x), and 

2^c^c 2(1) if y±>y°c 2(1) , 

> c )R C2(1) if y± < y° C2W , 



± 



(All) 



2(vr - c£±)i? f 



with <f>c = arcsin(x + /?? sin(a)) and y~ the y-coordinate of the center of the circle C2(l). The 

2(1) 



Lf=< 



2(pcRc 2{1) 



expression for Lf can be also written as 

if 0i G [@ ln {R,x);@t(R,x)}, 

2(n - <pt)Rc 2(1) if ©I € ]Q+(R,x);&g l (R,x)], 
where @f(R,x) is the value of the angle ©i at which x + = Rc 2 m- 

Considering R° < R < oo one can show that for x G [l;x°(R)} 



L s = < 



(A12) 






if 


0! 


G 


[0;© ln (i?,x)], 




Lf 


if 


0! 


G 


]e ln (R,x) ; e?™(R,x)], 


(A13) 


L b 


if 


0! 


G 


]e^(R,xy,e¥ x (R,x)}, 
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with @ ln (R, x) = a — 



(R,x), whereas for x € [x°(R);2r]] 



if 6iG [0;@ w (R,x)} 



L s = { 



if 9iG ]Q 10 (R,x);&?™(R,x)], 



(A14) 



if 0i £ ]e^ M (i2,x) ; e 



The probability s 2 (i) = s 2(i)(R, x ) defined by eq.(S2), can be obtained by integrating s'^JR, x, 0i) 



with respect to the angle ©i subject to the constraint that hb-arm 1 has already formed a bond (the 
distribution of possible orientations of hb-arm 1 assumed to be uniform) : 



Coefficient k^ 

The coefficient ks is found by assuming that hb-arms 1 and 2 of molecule S (see Figure 1) have 
already formed bonds in arbitrary orientations, determined by angles ©i and ©2 that they form with 
the x' axis and calculating the probability that the third hb-arm (say, hb-arm 3) can also form a 
bond subject to the constraint that the angle between any two hb-arms is equal to a. Certainly, for 
hb-arm 3 to form a bond its tip must not lie to the left of the plane LI (left boundary of the SHL). 
The orientations of hb-arms 1 and 2 are eventually averaged assuming their uniform distributions. 

Clearly, hb-arm 3 can form a bond with the same probability as in the bulk if the location of 
molecule S and orientation of its arms 1 and 2 (i.e., x, ©i, and ©2) are such that the tip of arm 3 is 
not located within the inner boundary of the particle SHL. Otherwise, hb-arm 3 cannot form a bond 




(A15) 



Therefore, according to eqs.(S3) and (S8), we have 




(A16) 
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at all. Therefore, s^ 2 ^ = s^ 2 ^ ®i) ©2), one can thus write 

63(2,1) if©i£ [eS^x^e^^x)] and e 2 € [eg;(i2,x,ei);eS(i2,x,ei)], 

S 3(2,l) ~~ * 

otherwise, 

(A17) 

(3) (3) 

where &\^(R,x) and @\J(R,x) are the minimum and maximum angles (both in the range from 
to 7r) that hb-arm 1 can have with the r-axis for three hydrogen bonds to form (for given R and x); 

(3) (3) 

Q 2x (R, x, ©i) and ©^ (R, x, @i) are the minimum and maximum angles (both in the range from 
to 7r) that hb-arm 2 can have with the t-axis for three hydrogen bonds to form simultaneously (for 
given R,x, and @S0R,x) < @i < ©SO^x) )• 

The mean probability s 3 / 2 ,i) that molecule 5 for given R and x forms a third hydrogen bond once 
its hb-arms 1 and 2 have already formed bonds can be obtained by averaging s^ 2 ^ over all possible 

orientations of its hb-arms 1 and 2, i.e., over ©i and ©2 (both angles distributed uniformly) 

etJ(R,x) e£>(fl,*,ei) 

/■ dgl /■ d©2 , 

S3{2 ' i} " y ©^(i?,x)-0 ln( i?,x) y e^w 7 ^^ 83 ' 2 ' 1 ' x ' 

(A18) 

Here Q ln (R,x) = H(x°(R) - x)Q ln (R, x) with H(x°(R)-x) the Heaviside step function; ©^(i^x) 

(2) 

and Q\J(R, x) are the minimum and maximum angles (both in the range from to tt) between hb- 

(2) 

arm 1 and r-axis with two hydrogen bonds formed simultaneously (for given R,x); 2 ; ; (i?,x,0i) and 

(2) 

2a . (-R, x, ©i) are the minimum and maximum angles (both in the range from to 7r) that hb-arm 
2 can have with the r-axis when two hydrogen bonds are formed simultaneously (for given R, x, and 
Q ln (R,x) < ©i < 0±l(R,x)). According to eqs.(S17) and (S3), the coefficient £3(2,1) {R, x) is thus 
given by 

k ( Rx )- [ & ^ {R ' X) dgi ^HR^e,) ^02 

3(2,1) ( > x > ~ J e( 3) {Rx) QM {R ^ x) _Q in{R ^ x) ] e ^ {Rxei) e 2x (R,x,e 1 )-@ 2n (R,x,@ l y 

(A19) 
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For R > R° and x G [0;x°(i?)], the angle Q$( R > 

x) is obtained as the solution of the equation 



(1 + cos{a))xl + (-1 + cos(a))y 2 ± + (-1 + cos(a))(z± -(R + rj + x)f = (A20) 

with respect to Oi (recall that x+, y±, and z± are all functions of R,x, and ©i). If < R < R° or 
x G ]x°(R),2rj\, the angle 

e^(R,x) = o. 



On the other hand, 



9 { 2(R,x) = Qf x (R,x) 



for any R and x. 

Next, for R > R° and x G [77; x°(R)], 

@f x (R,x) 

e%(R,x,e 1 ) = \ 



if e lG [e®(R,xy,e?f(R,x)], 



(A21) 



e£ ) (i?,x,e 1 ) = <( 



if Gi g [e 10 (i2,x) ; e^(fl,x)], 



(A22) 



7r — |7r — (a + 0i ) | otherwise, 
If < R < R° orxe ]x°(R),2r]], 

er x (R,x) 

-it — \tt — (a + ©i) | otherwise, 
Further, for R > R° and x G [77; x°(R)], 

§ - arcsin ( (^(i^)-^ if @1 € [eg(i2, x); 0^ M (i?,x)], 
|a — ©1 1 otherwise, 

where = Z2 n {R,x,@i) is the z-coordinate of the tip of hb-arm 3 when the tip of hb-arm 2 is 
located on the sphere S l g HL and Z s are the z-coordinate of the selected molecules S. 
For R > R° and x G [x°(R); x sep (R)], 

f- arcsin ( (^(^)-^) ) if @1 G [egjjfl, x); ej(fl, x)], 



©g(i?,x,0 1 ) = <( 



(A23) 



0g(i?,x,0 1 )= <( 



(A24) 



la — ©1 1 



otherwise, 
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where x scp = x scp (R) is the independent solution of a couple of simultaneous equations x + (R, x, 0i) = 

(3) (3) 

77sin(a/2) and y±(R,x,@\) = with respect to x, and @ llow (R,x) and (R, x) are the smaller and 
larger solutions of the equation x+(R, x, @±) = r/sin(a/2) (these solutions exist only at x < x sop and 

(3) (3) (3) 

at x = x scp they degenerate into a single solution @\ s ' ep {R) = @ llow (R, x sep (R)) = &\u P {R, x sep (R))). 
For any R > and x <G [x Bcp (R); 2rj], 

@^{R,x,@ 1 ) = \a-@ 1 \ if9i€ [0;Q? x (R,x)], (A25) 

For < R < R° and x € [r/; x scp (R)j, 

| - arcsin ( ^teg^l ) if 0i £ [egjjfl, x); 0^(12, x)], 



eg(i?,x,e 1 ) = < 



(A26) 



|a — ©i | otherwise, 



Coefficient 



Again, let us consider molecule S in the SHL of a spherical hydrophobe of radius R at a distance x 
from its surface. The coefficient k± is calculated by assuming that hb-arms 1,2, and 3 have already 
formed bonds and taking into account that, if the water molecule is far enough from the hydrophobic 
surface (but still in the LHS), hb-arm 4 can still form a bond if its tip is not within the inner boundary 
of the hydrophobe SHL; besides, the angle between any two of four hb-arms must be equal to a. 
While the orientations of arms 1 and 2 are arbitrary (with the angle between them equal to a), hb- 
arm 3 can have any of just two possible orientations determined by those of arms 1 and 2, whereas 
the orientation of hb-arm 4 is uniquely determined by the orientations of hb-arms 1, 2, and 3. (Again, 
the orientations of arms 1 and 2 are eventually averaged assuming their uniform distributions. 

What is the probability s' 4 , 3 2 ^ = s' 4 , 3 2 1 )(x / , Bi, 62) that hb-arm 4 will form a bond as well? 
Let us denote the angle between arm i {i = 1,2,3,4) and axis x' by Oj and introduce the same 
Cartesian coordinate system as above (in calculating k2(R,x)) but with the origin coinciding with 
S. The Cartesian coordinates of the tip of arm i (i = 1, 2, 3, 4) will be denoted by x,, yi, Z{. First of 
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all, in order for hb-arm 4 to form a bond, molecule S must be sufficiently far away from the inner 
boundary of the hydrophobe SHL (represented the sphere S l g HL ). The minimum distance x m i n , at 
which this is possible, depends on R and is equal to x°(R) defined above and determined by eq.(S9). 

Further, hb-arm 4 can form a bond with the same probability as in the bulk if the location of 
molecule S and orientation of its arms 1 and 2 (i.e., x',@\, and 62) are such that the tip of arm 4 
is not located within the sphere S l g HL . Otherwise, hb-arm 4 cannot form a bond at all. Keeping in 
mind that s^ 3 2 1) = s 4(3 2 x ; ®i> ®2), one can thus write 



S 4(3,2,l) 



63(2,1) if0 ie Af (4) ande 2 € Af (4) , 

(A27J 

otherwise. 

Here, Aw 4 , = A®, 4 , (R, x) is the range of angles (from to n) between hb-arm 1 and r-axis with 
four hydrogen bonds formed (for given R and x) and ^f(4) = ^2(4) x > ® x ) ^ s ^ ne ran S e °f angles 
(from to 7r) that hb-arm 2 can have with the r-axis when four hydrogen bonds form (assuming that 
01 € Af^(R,x) for given R,x). 

The mean probability 54(3 2,1) that a boundary molecule (i.e., molecule S) forms a fourth hydrogen 
bond once its hb-arms 1,2, and 3 have already formed bonds can then be obtained by averaging 
s 4(3 2 1) over au P oss ible orientations of its hb-arms 1 and 2, 0i and 62 (assumed to be distributed 
uniformly) : 

f ,10, f d02 , 

Sl( " 1,= a / eS^i-eSw,) / eg?Wx,e 1 )-eg(i i . I .e 1 ) S4,w,<B ' I ' ei,e2) - 

A l(4) A l(4) 

(A28) 

Thus, according to eqs.(S27) and (S3), the coefficient k±{R,x) is equal to zero for 77 < x < x°(R), 
whereas for x°(R) < x < 2rj it is given by 

fc 4 (3 2 / / -7^ • (A29) 

A l(4) A l(4) 
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A l(4) 



(A30) 



A 2(4) 



For any R > one can show that 

[0-MuL(R,x)] U [e ( S p (R,x);Qf x (R,x)} if x € [x°(R);x sep (R)}, 

[0;Q? x (R,x)] if xe [x sep (R);2ri}. 

For any R > and x € [x sep (.R); 277], 

[©S(©i);©^(©i)] if e lG [o ; e 10 (i2,x)] □ 
[0S(@i);@iil(i?,*)] U[0ii(^^);0if(^^)] if ©i e [eio^xjies^i)], 

(A31) 

where ©^(©i) = |a - ©i| and ©^©l) = tt - \ir - (a + ©i)| 
For any > and x <G [x°(.R); x sep (i?)], 

[0S(0i);02 4 J(0i)] if 0i e [o;e 10 (i?,x)] u [e^a^e^s)], 

if 01 G [0^(12,*); e®(fl,a;);], 

[02 4 n(0i);02 4 nl(^^ 5 0i)] U [6^,(i2, x, S^; 0^(i2, x)] otherwise 

(A32) 

(the last condition "otherwise" stands for @i € [0i O (i2, x); @u ow (R, %)] U [©itp(fl> <*0; 
The angles ©^(.R, x, 61)] and 2 ^ n (.R, x, 61)] in eqs.(S31) and (S32) are determined as 



^2(4) 



e$ x (R,x,e 1 )] = { 



arccos[z 2 (R, x, @i)] if |a - 0^ ep (i?)| = arccos[z 2 (R, x sop (R), of^R))], 



arccosfz^ (R, x, @i)] otherwise, 



(A33) 



and 



e£ ) J J R,x,e 1 )] = < 



arccos^ (i?, x, 61)] if 0^ ra (i?, x, 61)] = arccos[z 2 (i?, x, 0i)], 



(A34) 



arccos[z 2 x, @i)] otherwise. 
Here = z£(R, x, ©i) and = z 2 {R, x, @i) are two solutions of the quadratic equation 

%3\J 1 - (cos(a)/yi - (z 1 /y 1 )z 2 ) 2 - z\ = cos(a) - (cos(a)/yi - (z 1 /y 1 )z 2 )y3 - Z2Z3, (A35) 

with y\ = — sin 0i, z\ = cos©i, 2:3 = cos 0^(x, R), 2/3 = — (cos(a)/sin©i) + z% cot ©i, X3 = 

V 1 -y'i -4- 
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Captions 

to Figures 1 through 5 of the manuscript "Probabilistic approach to the length-scale de- 
pendence OF THE EFFECT OF WATER HYDROGEN BONDING ON HYDROPHOBIC HYDRATION" by 
Y. S. Djikaev and E. Ruckenstein. 

Figure 1. A water molecule in the surface hydration layer (SHL) of a spherical hydrophobic particle 
of radius R. The inner boundary of the SHL is a sphere of radius R + rj, while the outer (closer to 
the bulk water) boundary is shown as the sphere of radius R + 2r]. The molecule, shown as disk S, is 
located at the distance x from the hydrophobe surface. Two hb-arms of the molecule (arms 1 and 2) 
are in the plane of the Figure, while arms 3 and 4 are located out of the Figure plane (one of them 
under, the other above). The tips of hb-arms are shown as empty circles. The angle between any 
two hb-arms is a. The origin O of the Cartesian coordinate system lies coincides with the center of 
the spherical hydrophobe,. The angle between hb-arm i (i = 1,2,3,4) and axis x is denoted by Oj. 

Figure 2. The function n s (R, x) for a spherical hydrophobe in liquid water at temperature T = 293.15 
K (with rib = 3.65): a) n s vs £ = (xjr\ — 1) for various R's; b) n s vs R at different x's. 

Figure 3. The water density profiles near a hydrophobe of radius R (for R/r] = 1, 3, 5, 7, 10, 15, 30, 100) 
with e wp /e ww = 0.75 at T = 293.15 K and li = -11.5989 k B T. 

Figure 4. The water density profiles near a hydrophobe of radius R at T = 293.15 K and fj, = —11.5989 
ksT, with a) R/rj = 1 and b) R/r) = 1. Three curves each of Figs. 4a and 4b correspond to different 
£ wp /£„, with e wp /e ww = 0.6,0.75, and 1.2 for curves from bottom to top. 

Figure 5. The grand canonical free energy of hydration of a spherical hydrophobe of radius R. The 



30 



dimensionless AQ = AQ/(kBT(4irR 2 /r) 2 )) is plotted: a) vs R/rj for particles with e wp /e ww = 0. 
b) vs e wp /e ww for various i?'s (the curves are for R/rj = 1, 3, 7, 15,, and 30 from bottom to top) 
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